Astronomy & Astrophysics manuscript no. paper 


© ESO 2009 


June 7, 2009 





Electron density in the quiet solar coronal transition region from 
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Abstract 

Context. The sharp temperature and density gradients in the coronal transition region are a challenge for models and observations. 
Aims. We set out to get the average electron density (n c ) in the region emitting the S vi lines. We use two different techniques which 
allow to derive linearly-weighted (opacity method) and quadratically-weighted (Emission Measure method) electron density along 
the line-of-sight, in order to estimate a filling factor or to derive a thickness of the layer at the formation temperature of the lines. 
Methods. We analyze SoHO/SUMER spectroscopic observations of the S vi lines, using the center-to-limb variations of radiance, the 
center-to-limb ratios of radiance and line width, and the radiance ratio of the 93.3-94.4 nm doublet to derive the opacity. We also use 
the Emission Measure derived from radiance at disk center. 

Results. We get an opacity To at S vi 93.3 nm line center of the order of 0.05. The resulting average electron density (n e ), under simple 
assumptions concerning the emitting layer, is 2.4 • 10 16 irr 3 at T = 2 ■ 10 5 K. This value is higher than (and incompatible with) the 
values obtained from radiance measurements (2 • 10 15 irr 3 ). The last value leads to an electron pressure of 10~ 2 Pa. Conversely, taking 
a classical value for the density leads to a too high value of the thickness of the emitting layer. 

Conclusions. The pressure derived from the Emission Measure method compares well with previous determinations. It implies a 
low opacity of 5 10~ 3 to 10 -2 . The fact that a direct derivation leads to a much higher opacity remains unexplained, despite tentative 
modeling of observational biases. Further measurements, in S vi and other lines emitted at a similar temperature, need to be done, and 
more realistic models of the transition region need to be used. 

Key words. Sun : atmosphere - Sun: transition region - Sun: UV radiation 



1. Introduction 

In the simplest description of the solar atmosphere, where it is 
considered as a series of concentric spherical layers of plasma at 
different densities and temperatures, the transition region (here- 
after TR) between the chromosphere and the corona is the thin 
interface between the high-density and low-temperature chro- 
mosphere (a few 10 16 irr 3 hydrogen density at about 10 4 K) and 
the low-density and high-temperature corona (about 10 14 irr 3 
at 10 6 K). The variation of temperature T and electron number 
density n e has been mostly derived from the modelling of this 
transition region , where radiative losses are balanced by thermal 
conduction (e.g. lMariskalll993l:lAvrett & Loeserll2008l) . 

Measurements of the electron density usually rely either on 
estimation of the Emission Measure or on line ratios. On one 
hand, using absolute line radiances, the Emission Measure (EM) 
and Differential Emission Measure (DEM) techniques provide 
(/ig) at the formation temperature of a line (or as a function of 
temperature if several lines covering some range of temperatures 
are measured). On the other hand, the tec hnique of lin e radiance 
ratios provides a wealth of values of n e dMasonll 19981) with the 
assumption of uniform density along the line-of-sight, and with 
an accuracy limited by the accuracy of the two respective ra- 
diance measurements: typically, a 15 % uncertainty on line ra- 
diance measurement leads to 30 % uncertainty on the line ratio 
and then to about a factor 3 uncertainty on the density. However, 
for a given pair of lines, this technique only works in a limited 
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range of densities. Let us add that the accuracy is also limited by 
the precision of atomic physics data. 

Here we propose to use also the concept of opacity (or op- 
tical thickness) in order to derive the population of the low (ac- 
tually the ground) level i of a given transition i — > j, and then 
the electron density. At a given wavelength, the opacity of a col- 
umn of plasma corresponds indeed to the sum of the absorp- 
tion coefficients of photons by the individual ions in the column. 
The op acity can be derive d by different complementary tech- 
niques dDumont et alj 19831) if many measurements are available 
with spatial (preferably center-to-limb) information. This is the 
case in a full-Sun o bservations program by the SoHO /SUMER 
UV s pectro-imager dWilhelm et alj] 1 995t iPeterl 1 1 999h Peter & 



Judge ll999l) run in 1996. In particular, thanks to a specific "com- 
pressed" mode, a unique dataset of 36 full-Sun observations in 
S vi lines has been obtained; this makes possible to derive at the 
same time (« e ) from opacity measurements and {n\) from line 
radiance measurements (via the EM). 

We have already used this data set in order to get properties 
of turbulen ce in t he TR dBuchlin et al.ll2006l). Note that here, 



contra ry to IPeterl dl999h : |Peter & Judge! (119991) : iBuchlin et all 



(120061) . we are not interested in the resolved directed velocities 
or in the non-thermal velocities but in the line radiances, peak 
spectral radiances and widths. Also note t hat, along with the 
modelling work of lAvrett & Loeserl d2008l) . we do not distin- 
guish network and internetwork (anyway a difficult task at the 
limb) and aim at a precise determination of the properties of an 
average TR. 

This paper is organized as follows: we first present the 
data set we use, then we determine opacities and radiances of 
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Figure 1. Raw line profiles from the context spectrum taken on 

4 May 1996 at 07:32 UT at disk center with an exposure time of 
300 s. The profiles are averaged over pixels 50 to 299 along the 
slit (1 x 300 arcsec, detector A), with no prior destretching of the 
data. 

5 vi 93.3 nm, we get two determinations of density in the region 
emitting the S vi 93.3 nm line, we discuss the disagreement be- 
tween the two determinations (especially possible biases), and 
we conclude. 



2. Data 

2.1. Data sets 

We use the data from a SoHO/SUMER full-Sun observation pro- 
gram in S vi 93.3 nm, S vi 94.4 nm and Ly s designed by Philippe 
Lemaire. The spectra, obtained with detector A of SUMER and 
an exposure time of 3 s, were not sent to the ground (except for 
context spectra) but 5 parameters ("moments") of 3 lines were 
computed on-board for each position on the Sun: 

- (1) peak spectral radiance, (2) Doppler shift, and (3) width 
of the line S vi 93.3 nm, 

- (4) line radiance (integrated spectral radiance) of the line 
Lye 93.8 nm, 

- (5) line radiance of the line S vi 94.4 nm. It must be noted 
that this line is likely to be blended with Si vm. 

The detailed characteristics of these lines can be found in 
Table Q] A list of the 36 observations of this program run 
throughou t year 1996, close to s olar minimum, can be found in 
Table 1 of iBuchlin et~aT1 d2006b . These original data constitute 
the main data set we use in this paper, hereafter DS 1 . They are 
complemented by a set of 22 context observations from the same 
observation program, that we use when we need the full profiles 
of the spectral lines close to disk center: the full SUMER de- 
tector (1024 x 360 pixels) has been recorded at a given position 
on the Sun at less than 40 arcsec from disk center and with an 
exposure time of 300 s. This data is calibrated using the Solar 
Software procedure sum_read_corr_fits (including correc- 
tion of the flat field, as measured on 23 September 1996, and of 
distortion), and it will hereafter be referred to as DS2. 

2.2. Averages of the data as a function of distance to disk 
center 

In order to obtain averages of the radiances in data set DS 1 as a 
function of the radial distance r to the disk center, and as a func- 
tion of /i, the cosine of the angle between the normal to the solar 
"surface" and the line-of-sight, we apply the following method, 
assuming that the Sun is spherical: 



- We detect the limb automatically by finding the maximum 
of the S vi 93.3 nm radiance at each solar-y position in two 
detection windows in the solar- x direction, corresponding to 
the approximate expected position of the limb. This means 
that the limb is found in a TR line and is actually approxi- 
mately 3 arcsec above the photosphere. However, this is the 
relevant limb position for the geometry of the S vi 93.3 nm 
emission region. 

- We fit these limb positions to arcs of a circle described by 
x(y) functions, and we get the real position (a, b) of the solar 
disk center in solar coordinates (x,y) given by SUMER, and 
the solar radius Rq (this changes as a function of the time of 
year due to the eccentricity of SoHO's orbit around the Sun). 
The solar radius is evaluated for the observed wavelength of 
93.3 nm. 

- We choose to exclude zones corresponding to active regions, 
as the aim of this paper is to obtain properties of the TR in 
the Quiet Sun. 

- For each of the remaining pixels, we get values of the radial 
distance r = y(x — a) 2 + (y - b) 2 to disk center and of /v = 
VI - (r/R Q ) 2 . 

- We compute the averages of each moment (radiances and 
widths) in bins of r/RQ and in bins of 1 /fi. 

The resulting averages as a function of r/RQ and of 1 //i are plot- 
ted in Fig. |2] (except for the S vi 93.3 nm Doppler shift, which 
will not be used in this paper). The radiances are approximately 
linear functions of l/fi for small 1 /fi, as expected from optically 
thin lines in a plane-parallel geometry. Such a behavior actually 
validates the consideration of a "mean" plane-parallel transition 
region, at least for l/fx < 10 or 6 < 84°. 

3. Determination of opacities 

3.1. Using center-to-limb variations 

We follow here the method A proposed bv lDumont et alj dl983l) . 
Assuming that the TR is spherically symmetric and that it can be 
considered as plane-parallel when not seen too close to the limb, 
that the lines are optically thin, and that the source function S 
is constant in the region where the line is formecQ, the spectral 
radiance is: 

/ Gu) = 5(l-exp(-T //i)) (1) 

where the subscript is for the line center and r is the opacity of 
the emitting layer at disk center. Then: 



I (M) = 7 (1) 



1 - exp(-T //j) 
1 - exp(-ro) 



(2) 



and a fit of the observed /oGu) by this function, with Io(l) and to 
as parameter^, gives an estimate of to. 

For the lines for which only the line radiance E is known 
(S vi 94.4 nm and Ly s), we need to fit this function, with t and 
£(1) as parameters^]: 



E(M) =E(l) 



l(l-exp(-K" ; )) du 
j;(l-exp(-T 0e -" 2 )) du 



(3) 



1 We release this assumption in Sec.[5] 

2 Note that, contrary to Dumont et afj dl983h . we take /o(l) as an ad- 
ditional parameter. This is because by doing so, we avoid the sensitivity 
of / (1) to structures close to disk center, and because the first data bin 
starts at 1/yu = 1 instead of being centered on 1/fi = 1 

3 We take here £(1) as a parameter for the same reason as we did 
before for /o(l). 
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Table 1. Spectral lines present in the data sets, with parameters computed by CHIANTI and given by previous observations. 



Ion Transition j — > i 


CHIANTI" 

log r max (K) Wavelength (A) Radiance^ 


Curdt et al. (2001) 
Wavelength (A) Peak radiance* 


Svi 2p b 3p 2 P 3/2 -» 2p b 3s 2 Si /2 
Lye 6p 2 P 3/2 -> ls 2 Si /2 
Si vm 2s 2 2p 3 2 P 3/2 -> 2s 2 2p 3 4 S 3/2 
Svi 2p 6 3p 2 P, /2 -> 2p 6 3s 2 S 1/2 


5.3 933.3800 3.81 ■ 10" J 

5.9 944.4670 5.24 • 10~ 3 
5.3 944.5240 1.91 ■ 10" 3 


933.40 0.57 
937.80 1.07 
944.34 0.14 
944.55 0.29 



" Using the "Arnaud & Raymond" ionization fractions file, the "Sun coronal" abundance file and the "Quiet Sun" DEM file. CHIANTI does not 
include data for the Hydrogen lines (Ly e in particular). 

Radiances are given in Wm sr , and peak spectral radiances are given in Wm^sr'nm -1 . 
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Figure 2. Average of the data as a function of r//?© (top panels) 
and as a function of l//i (bottom panels). 



These theoretical functions of fi are then plotted for differ- 
ent values of the parameter to over the observations in Fig. [3] 
for all three lines (either for the peak spectral radiance or the 
line radiance, depending on the data). We have performed a non- 
linear least-squares fit using the Levenberg-Marquardt algorithm 
as implemented in the Interactive Data Language (IDL); it gives 
the parameter To. The uncertainties on each point of the E(p) or 
functions (an average on Nd pixels) that we take as input 
to the fitting procedure come mainly from the possible presence 
of coherent structures such as bright points: the number of such 
possible structures is of order Nd/N s , where N s is the size of 
such a structures (we take N s = 100 pixels), and then the uncer- 
tainty on / or E is crj y/Nd/N s where <x is the standard deviation 
of the data points (in each pixel of a l//u bin). Compared to this 
uncertainty, the photon noise is negligible. 

The results of the fits on the interval l/fj. e [1,5] are shown 
in Fig. [3] as far as to is concerned, they are 0.1 13 for moment (1) 
(S vi 93.3 nm peak spectral radiance) and 0.244 for moment (5) 
(S vi 94.4 nm radiance, blended with Si vm). The approximations 
we used in writing Eq. (fl} are not valid for the optically thick 
Ly s line, hence the bad fit. On the other hand, these approxi- 
mations are valid for both the S vi lines, as long as 1 /// is small 
enough. For large 1//J. there is an additional uncertainty resulting 
from the determination of the limb. 

These results are somewhat sensitive to the limb fitting: a 
10~ 3 relative error in the determination of the solar radius leads 
to a 7 10~ 2 relative error on To. As 10~ 3 is a conservative upper 
limit of the error on the radius from the limb fitting, we can con- 
sider that 7 10~ 2 is a conservative estimate of the relative error 
on To resulting from the limb fitting. 



3.2. Using center-to-limb ratios of Svi 93.3 nm width and 
radiance 

The variation with position of the Svi 93.3 nm line width 
(see Fig. |2]l can be interpreted as an opacity sat uration of the 
Svi 93 .3 nm line at the limb, and then method B of Dumont et al. 
(1983) can be applied. This method relies on the measurement 
of the ratio d = AA*/AA* C of the FWHM at the limb and at the 
disk center: the optical thickness at line center ?o at the limb is 
given by solving 



This expression comes from Du mont et al.l d 19831) and assumes 
a Doppler absorption profile exp(-M 2 ) with u = AA/AAo- Here, 
contrary to the case of the peak spectral radiance ratio, the func- 
tion and its derivative with respect to to and £(1) cannot be 
computed analytically anymore, and we need to estimate them 
numerically; this is done by a fast method, using a Taylor ex- 
pansion of the outermost exponential of both the numerator and 
denominator of Eq. (|3). 



2(l -exp(-f e- rf2ln2 )) = 1 -exp(-fo) 



(4) 



(this is Eq. 4 of iDumont et alj[l983l where a sign error has been 
corrected) and then the opacity at line center to is given by solv- 
ing 

7o0"=l) l-exp(-To) 



I Q (/u = 0) 1 - exp(-fo) 



(5) 
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Figure 3. Diamonds: average profiles of the radiance data (mo- 
ments 1, 4 and 5) as a function of 1 /fi, normalized to their values 
at disk center. Dotted lines: theoretical profiles for different val- 
ues of to. Plain lines: fits of the theoretical profiles to the data, 
giving the values for t : 0.1 13 for (1) and 0.244 for (5). The fit 
for Ly e is bad because this line is optically thick. 

Using the full-Sun S vi 93.3 nm compressed data set DSlQ, 
we find that the ratio d is 1.274 and then to is 1.53. Finally, 
we use the Svi 93.3 nm peak spectral radiance ratio Iq(ji = 
1)1 Mil = 0) = 0.062 to get r = 0.05. 

4 Although not obvious from the data headers, moment (3) corre- 
sponds to the deconvoluted FWHM of Svi 93.3 nm, as is confirmed 
by a comparison with the width obtained from the full profiles in 
data set DS2 and deconvoluted using the Solar Software procedure 
con_width_4. 



3.3. Using the S vi 94.4 - 93.3 line ratio 

The theoretical dependence of the S vi 94.4 - 93.3 peak radiance 
line ratio as a function of the line opacities and source functions 
is: 

„ A),933 S933 1 - exp(-To,933) ,„ 

K — — (o) 

2 7o,944 2 ^944 1 - exp(-To,944) 

For this doublet, we assume S 933 = S 944 and ^933 = 2^944 (be- 
cause the oscillator strengths are in the proportion ^33 = 2/944). 
Then K reduces to 

1 1 - exp(-T ,933) 1 + exp(-T ,933 /2) ,_. 

A = — ( / ) 

2 1 -exp(-T 0; 933 /2) 2 
and we get to,933 from the observed value of K: 

t-0,933 = -2 ln(2tf-l) (8) 

The difficulty comes from the S vi 94.4 nm blend with the 
Si viii line. In order to remove this blend, we have analyzed the 
line profiles available in data set DS2. After averaging the line 
profiles over the 60 central pixels along the slit, we have fitted 
the S vi 93.3 nm line by a Gaussian with uniform background and 
the S vi 94.4 nm line blend by two Gaussians with uniform back- 
ground. We have then computed the Gaussian amplitude from 
these fits for both S vi lines, and this gives 7o,933 and 7o,944, and 
then K, that we average over all observations. From this method 
we get to,933 = 0.089. 

The same kind of method could in theory be used for the S vi 
94.4 - 93.3 line radiance ratio 

K _ £933 _ S933 i(l-exp(-T ,933 g-" 2 ))d M 
2 £944 2 S944 j R (l - exp (-t 0; 944 e~" 2 )) du 

with, again, S 933 = 5 944 and 7-0,933 = 2to,944. As for method A, 
the integral makes it necessary to invert this function of ^933 
numerically, in order to recover 70,933 f° r a given observed value 
of K. As K is decreasing as a function of To,933, this is possible 
by a simple dichotomy. However, the average K from the obser- 
vations is greater than 1, which makes it impossible to invert the 
function and get a value for To. 

3.4. Discussion on opacity determination 

It is clear that the three methods provide different values of 
the opacity at disk center. We confirm the result of Dumont 
et al. d!983h . obtained in different lines, by which the method 
of center- to-li mb ratios of width and radiance (Sec. 13.21 or 
method B in Dum ont et al.l fl983) provides the smallest value 
of the opacity. As mentioned by these authors, the center-to- 
limb variations method (Sec. 13. II or meth od A) overestimates th e 
opacity for different reasons described in Dumo nFet al.l JT983), 
among which the curvature of the layers close to the limb and 
their roughness. The method of line ratios (Sec |3.31 or method C) 
also provides larger val ues of the opaci t y, alth ough free from ge- 
ometrical assumptions; Dum ont et alj d!983h interpret them as 
resulting from a difference between the source functions of the 
lines of the doublet. 

This does not mean that there are no additional biases. For 
instance, we have adopted a constant Doppler width from center 
to limb; actually this is not correct since at the limb the observed 
layer is at higher altitude, where the temperature and turbulence 
are higher than in the emitting layer as viewed at disk center. 
Consequently, the excessive line width is wrongly interpreted 
as only an opacity effect. However, it seems improbable that a 
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27.4% increase of Doppler width from center to limb can be en- 
tirely interpreted in terms of temperature (because of the square- 
root temperature variation of Doppler width) and turbulence (as 
the emitting layer is — a posteriori — optically not very thick). 



4. First estimates of densities 

4.1. Densities using the opacities 

The line-of-sight opacity at line center of the S vi 93.3 nm line is 
given by 



used to estimate (n 2 .) As in the emitting region of thickness As, 



TO 



I 



k Vo nsm,i( s ) ds 



(10) 



where the integration is along the line-of-sight. The variable 
«s vi,; is the numerical density of S vi in its level i, which can 
be written as 



"Svi,; «Svi ,,„,"H 

«Svii = Abund(S) — n e 

«Svi «S «e 



(11) 



where Abund(S) = «s/«h is the Sulfur abund ance in the corona 
(IP" 4 - 73 according to the CHIANTI database. iDere et alj[l997t 
lLandi et ai1l2006l) «svi,;/"Svi is the proportion of S vi at level i, 
ttSvi/ns is the ionization fraction (known as a function of temper- 
ature) and rtn/n e - 0.83 is constant in a fully ionized medium 
as the upper transition region. In this work i is the ground state 
i = 1, and as in this region «svi,i/«Svi is very close to 1, we 
will drop this term from now. The variable lc Vo is the absorption 
coefficient at line center frequency vq for each S vi ion, given by: 



hvn 
k - -R. 



1 



47T -y/jrAvo 



(12) 



where B,j is the Einstein absorption coefficient for the transition 
i -> j (i.e., 2p 6 3s 2 Si /2 -> 2p 6 3p 2 P 3/2 ) at A = 93.3 nm and 
integration over a Gaussian Doppler shift distribution has been 
done (Avo is the Doppler width in frequency). Using: 



d _ 8j _ 8j Ajj 
Bij - —Bji 



gi gilhv^/c 2 
with gj/gi = 2 and Aq = c/vo, this gives: 



4tt 3 / 2 c AA d 



(13) 



(14) 



Finally, for an emitting layer of thickness As and average 
electron density (n e ), we have: 

tfAfi "Svi ,,„ «h, , K > 

To = am? a i Abund(S)— <n e > As (15) 

Taking r = 0.05, we get <n e )As = 4.9 ■ 10 21 nr 2 . Then, with 
As = 206 km (the altitude interval corresponding to the FWHM 
of the S vi 93.3 nm contribution function G(T) as computed by 
CHIANTI), this gives <n e ) = 2.4 ■ 10 16 nr 3 . 

4.2. Squared densities using the contribution function 

The average S vi 93.3 nm line radiance at disk center obtained 
from data set DS2 (excluding the 5% higher values which are 
considered not to be part of the quiet Sun) is £ = 1.4 ■ 
10 _2 Wrrr 2 sr _1 (to be compared to the value 3.81 ■ 10~ 3 given by 
CHIANTI with a Quiet Sun DEM — see Table Q). This can be 



E = J GiT{s)) » e 2 (s) ds * G((T}) <« 2 > As (16) 

where G{T) is the contribution function and the integral is on 
the line-of-sight and where we have made the assumption that 
to <k 1 . We take the average temperature in the emitting region 
to be (T) = r raax = 10 5 ' 3 K, and, for densities of the order of 
10 16 nr 3 , the gofnt function of CHIANTI gives G«r» = 1.8 • 
10~ 37 Wm 3 sr 1 . We finally get 



<« 2 >As = 8.4- 10 35 nT 5 



(17) 



With again As = 206 km, we get (n e )RMS = 2.0 • 10 15 nr 3 . 
Assuming an uncertainty of 20% on E, the uncertainty on 
(«e>RMS would be 10% for a given As. 

5. Discussion of biases in the method 

One of our aims when starting this work was to determine a fill- 
ing factoi0 

<«e) 2 



<» 2 > 



(18) 



in the S vi-emitting region. This initial objective needs to be 
revised, since we get / = 144, an impossible value as it is 
more than 1. Our values of d ensities can be compared to the 
density at log T = 5.3 in the lAvrett & Loeseri d2008l) model 
( 1 .7 • 10 15 nr 3 ): our value of <n e ) is an order of magnitude higher, 
while (n e )RMS = V^"e) is about the same (while it should be 
higher than (n e ))- Our value of intensit y is compatible with av- 
erage values from other sources, such as De l Zanna et al.l d2001h 
(see their Fig. 1). 

Given the same measurements of to and E, one can instead 
start from the assumption of a filling factor / e [0, 1] and deduce 
As: 



1 «»e)As) 2 

As = 

/ <» 2 >As 



(19) 



where the numerator and denominator of the second fraction are 
deduced from Eq. d 1 5b and ( 116b respectively. With the values 
from Sec. |4] this gives As > 29 Mm = Q.Q4Rq, a value much 
larger than expected. 

In any case, there seems to be some inconsistencies around 
log T — 5.3 between our new observations of opacities on one 
hand, and transition region models and observations of intensi- 
ties on the other hand. We propose now to discuss the possible 
sources of these discrepancies, while releasing, when needed, 
some of the simplistic assumptions we have made until now. 

5. 1 . Assumption of a uniform emitting layer 
5.1 .1 . Bias due to this assumption 

When computing the average densities from the S vi 93.3 nm 
opacity and radiance, we have assumed a uniform emitting layer 
at the temperature of maximum emission and of thickness As 
given by the width of contribution function G(T). However, 
the different dependences in the electron density of Eqs. d 1 0b 
and (fTBT l — the first is linear while the second is quadratic — 
means that the slope of the n e (s) function affects differently the 
weights on the integrals of Eqs. d 1 0b and (fl6]i: a bias, different for 



5 We explain this definition of the filling factor in Appendix lAl 
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To and E, ca n be expected, and here we explore this effect start- 
ing from the lAvrett & Loeserl d2008h model, which has the merit 
of giving average profiles of temperature and density (among 
other variables) as a function of altitude s. 



Opacity. Using the lAvrett & Loeser! d2008) profiles and atomic 
physics data, we get to = 0.008. Then, using the same simplistic 
method as for observations (still with a uniform layer of thick- 
ness As - 206 km), we obtain (« e ) = 2.4 ■ 10 15 irT 3 , a value 
only 40% higher than the density at logT = 5.3 in this model 
(1.7 ■ 10 15 m" 3 ). 



Radiance. Using the same lAvrett & Loeserl d2008l) profiles and 
the CHIANTI contribution function G(T), we get £ = 1.3 • 
10- 2 Wm- 2 sr _1 . 

Then, using the same simplistic method as for observations, 
we obtain (n e )RMS = 1.9 ■ 10 15 irT 3 , a value 12% higher than the 
density at log T = 5.3 in this model. 

We see then that the assumption of a uniform emitting layer 
has a bias towards high densities, which is stronger for the opac- 
ity method than for the radiance method. A filling factor com- 
puted from these values would be / = 1.5, while it has been 
assumed to be 1 when computing m and E from the Avrett & 
Loeser (2008) model: this can be one of the reasons contributing 
to our too high observed filling factor. 

This differential bias acts in a surprising way as, due to the ra 2 
term in Eq. ([TBI one would rather expect the bias to be stronger 
for E than for to; however, it can be understood by comparing 
the effective temperatures for tq and E, which are respectively: 



'eff,r - 



T e ffE - 



1Q 5.38 K 



/ T(s)n e (s)K(T(s)) ds 

J n e (s)K(T(s)) ds 
jT(s)nl(s)G(T(s)) ds_ = ^ r 
Jnl(s)G(T(s)) ds 



(20) 
(21) 



where K(T) - k Vn (T) ns n /n e , while T(s) and n e (s) are from 
I Avrett & Loeserl (|2008). The higher effective temperature for E 
than for to means that the bias is more affected by the respective 
shapes of the high-temperature wings of G(T) and K(T) than by 
the exponent of n e in the integrals of Eqns. ( TTZl ) and dl6t . 

It can be pointed out here that the difference between the 
K(T) and G(T) kernels lies in the fact that G(T) (unlike K(T)) 
not only takes into account the ionization equilibrium of S vi, but 
also the collisions from i to j levels of S vi ions. 



5.1 .2. Releasing this assumption: a tentative estimate of the 
density gradient around log T = 5.3 

In Sec. |5.1| we have incidentally sho wn that the radiance com- 
puted with the lAvrett & Loeserl d2008l) profiles and the CHIANTI 
contribution function G(T) is a factor 3 higher than the radi- 
ance computed directly by CHIANTI using the standard Quiet 
Sun DEM (see Table Q]). This is simply because the DEM com- 
puted from the temperature and density profiles of the Avrett & 
Loeser d2008h model is different than the CHIANTI DEM, as 
can be seen in Fig. [4] In particular, the lAvrett & Loeserl f2008) 
DEM is missing the dip around log T = 5.5 that is obtained from 



6 The reason for this is that the Avrett & Loesei| d2008l) model is 
determined from theoretical energy balance and needs further improve- 
ments in order to reproduce the observed DEM (E. Avrett, private com- 
munication). 



most observations; at log T = 5.3 it is a factor 3 higher than the 
CHIANTI Quiet Sun DEM. 

We model the upper transition region locally around log To = 
5.3 and so = 2.346 M m (chosen because T(so) = Tq in the 
lAvrett & Loeserl I200 8 model) by a vertically stratified plasma 
at pressure Po = 1.91/io&b7o (we consider a fully ionized coro- 
nal plasma) and: 



T(s) 



no 
n e (s) 



s — st 
s - s T 



for s > st 



(22) 



These equations were chosen to provide a good approximation 
of a transition region, with some symmetry between the oppo- 
site curvatures of the variations of T and n e with altitude. The 
parameters of this model atmosphere are the pressure Po and st 
(with st < so), which can be interpreted as the altitude of the 
base of the transition region. Given the constraint T{so) = To 
that we imposed when building the model, with To and so fixed, 
st actually controls the derivative of T(s) at s — so'- 



T'(s ) 



To 



2{s - s T ) 



or st 



■so 



T 



2T'(s ) 



(23) 



We plot in Fig. [5] some temperature profiles from this simple 
transition region model, for different model parameters T'(sp) 
(Po on ly affects the scale of n e (s)). For the lAvrett & Loeserl 
(2008) model, P„ = 8.7 • 10~ 3 Pa and T'(s ) = 0.45 KirT 1 , and 
the corresponding model profile is also shown. 

We propose to use such models along with atomic physics 
data and the equations of Sec.|4]to compute To and £ as a func- 
tion of model parameters Po and T'(so), as shown in Fig. [6] As 
the slopes of the level lines are different in the to(Po, T'(sq)) 
and E(Po, T'(sq)) plots, one would in theory be able to estimate 
the parameters (Po, T'(so)) of the best model for the observation 
of (To.obs, Eobs) by simply finding the crossing between the level 
lines t (Po, T'(s )) = T ,ob s and E(P , T'(s )) = E ohti - 

In practice however, the level lines for our observations of 
to and E do not intersect in the range of parameters plotted in 
Fig- El corresponding to realistic values of the parameters. As a 
consequence, it is not possible to tell from these measurements 
(from a single spectral line, here S vi 93.3 nm), what is the tem- 
perature slope and the density of the TR around the formation of 
this line. 

If we now extend the range of T'(sq) to unrealistically 
low values, a crossing of the level lines can be found below 
logP = -3.5 and T'(s ) = 5 mK/m. Given the width of G(T) 
for S vi 93.3 nm, this corresponds to As > 20 Mm, a value con- 
sistent with the one obtained from Eq. d 1 9t and which is also 
much larger than expe cted. 

Let us note that iKeenanl d!988l) derived a much lower 
S vi 93.3 nm opacity value (to = 1.1 10 -4 at disk center ) from a 
comp utation implying the cells of the network model of Gabriel 
(119761) . However, while our value of to seems to be too high, the 
level lines in Fig. [6] show that an opacity value to = 1.1 10~ 4 
would be too low: from this figure we expect that a value com- 
patible with radiance measurements and with realistic values of 
the temperature gradient would be in the range 5 10~ 3 to 10 -2 . 

5.2. Anomalous behavior of Na-like ions 
Followin g works such aslDupred dl972 ) for Li-like ions. Judge 



et al. (1 19951) report that standard DEM analysis fails for ions 
of the Li and Na is oelectronic sequences; in particular, for S vi 
(which is Na-like), iDel Zanna et al.l d200lh find that the atomic 
physics models underestimate the Svi 93.3 nm line radiance E 
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Figure 4. Quiet Sun standard DEM from CHIANTI (p lain line) 
and DEM computed from the lAvrett & Loeserl d2008h temper- 
ature and density profiles. The dotted lines give the DEMs for 
log T = 5.3, the maximum emission temperature of the S vi 
lines. 
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Figure 5. Temperature as a function of altitude in our local 
transition region simple models aroun d Tp = 10 53 and sp = 
2.346 Mm. The temperature profile from lAvrett & Loeserl f2008) 
is shown with the diamonds signs, and the simple model with 
the same temperature slope is shown with a dashed line. 

by a factor 3. This fully explains the difference between our ob- 
servation of E and the value computed by CHIANTI (Table [TJ. 
However, this means also that where G(T) from CHIANTI is 
used, as in Eq. (TToT l. it presumably needs to be multiplied by 3. 
As a result, one can expect (n e )RMS to be lower by a factor 1.7, 
resulting into a filling factor of 415 (actually worse than our ini- 
tial result). 

The reasons for the anomalous behavior of these ions for 
G(T), which could be linked to the ionization equilibrium or to 
collisions, are still unknown. As a result, it is impossible to tell 
whether these reasons also produce an anomalous behavior of 
these ions for K{T), hence on our measurements of opacities 
and on our estimations of densities: this could again reduce the 
filling factor. 

5.3. Cell-and-network pattern 

When analyzing our observations, we have not made the dis- 
tinction between the network lanes and the cells of the chro- 
mospheric supergranulation. Here we try to evaluate the effect 
of the supergranular pattern on our measurements, by using a 




-3.5 -3.0 -2.5 -2.0 -1.5 -1.0 

log P [Pa] 



0.8 - 



E 0.6 - 



0.2 - 




-3.5 -3.0 -2.5 -2.0 -1.5 -1.0 

log P [Pa] 

Figure 6. S vi 93.3 nm opacity to (top panel) and line radiance E 
(middle panel) as a function of model parameters Po and T'(sp). 
The level lines close to our actual observations are shown as 
plain lines for To and as dashed lines for E. The bottom panel 
reproduces these level lines together in the same plot. The pa- 
rameters (Po, T'(sq)) estimated from the lAvrett & Loeserl d2008) 
model at T — Tp are shown with the diamond sign on each plot. 



2D model emitting layer with a simple "paddle wheel" cell-and- 
network pattern: in polar coordinates (r, 0), the emitting layer is 
defined by Ri < r < R2; in the emitting layer, the network lanes 
are defined by 9 e [0, 68] mod AO and the cells are the other 
parts of the emitting layer, with A6 the pattern angular cell size 
(an integer fraction of 2n) and 59 the network lane angular size. 
The network lanes and cells are characterized by different (but 
uniform) source functions S , densities «s V i and absorption coef- 
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of Pr(or): 



rr cos ads 

(t (8, Pr (a))>„ = n e (s)K(T(s))—- Pr(a) da (25) 

J J cos(# + a) 

= Io(£2!^\ ,^ mMa)) (26) 
fx \cos(8 + a)l a /j 

The opacity to = To/jU is corrected by the factor /?(#, Pr(or)) de- 
fined in the previous equation. We recover ft — 1 for Pr(a-) = 
6(a), i.e., when there is no roughness. 

We immediately see that f3 = 1 for 9-0, for any Pr(o'): 
roughness (as modelled here by incompressible vertical dis- 
placements) does not change the optical thickness at disk cen- 
ter. Nevertheless, the estimate of optical thi ckness at disk center 
from observations in Sec. l3.1l (method A of Du mont et al.l fl983) 
is affected by this roughness effect. 

Coming back to (to), we take Pr(a) = cos 2 (7ra/2A)/A, and 
we compute /3 numerically (A represents the width of Pr(a) and 
can be thought as a quantitative measurement of the roughness). 
The results, shown in Fig.|9j indicate for example that the mod- 
elled roughness with A = n/5 increases the opacity by 9% at 
= 1.5 (corresponding to 9 = 45°). This is a significant ef- 
fect, and we can evaluate its influence on the estimate of to in 
Sec. 13. II in the theoretical profiles of Io(p) and E(fj) (Eq. |2]-[3]), 
To/fi needs to be replaced by Toln X /3. As /3 > 1 for a rough 
corona, this means that the value of to determined from the fit 
of observed radiances to Eq. ©-([3]) is overestimated by a fac- 
tor corresponding approximately to the mean value of ft on the 
fitting range. 

In this way, we have given a quantitative value for the over- 
estimation factor of to by the method of Sec. 13. II thus extending 
the qu alitative discussion on roughness found in iDumont et al.l 
(Il983h . This factor, of the order of 1.1 may seem modest, but 
one needs to remember that the fit for obtaining to in Sec. 13.11 
was done on a wide range (l/fi from 1 to 5, or 9 from to 78 
degrees) that our roughness model cannot reproduce entirei)Q. 

One can think of different roughness models representing the 
strong inhomogeneity of the TR, f or instance w i th a d ifferent 
and very peculiar roughness model iPecker et alj {1988) obtain 
an overestimation factor of more than 10 under some conditions. 
This means that our values of to may need to be decreased by a 
large factor due to a roughness effect. 

Roughness models can be seen as simplified models of the 
fine structure of the TR, which is known to be heterogeneous 
at small scales. Indeed, in addition to the chromospheric net- 
work pattern that we have already modelled in Sec. 15.31 the 
TR contains parts of different structures, with different plasma 
properties, lik e the base of large l oops and c oronal funnels, 
smaller loops dDowdv et alj 1 1986t lPeterll2001h . and spicules. 
Furthermore, the loops themselves are likel y to be composed of 
strands, which can be hea ted independently (Cargill & Klimchuk 
20Q4; |Parenti et al] |2006). The magnetic field in these structures 
inhibits perpendicular transport, and as a consequence the hori- 
zontal inhomogeneities are not smoothed out efficiently. 

6. Conclusion 

We have first derived the average electron density in the TR from 
the opacity to of the S vi 93.3 nm line, obtained by three differ- 



Figure 7. Average spectral radiance at line center Io as a function 
of \/n for a uniform layer (dashed line) and for a model layer 
with a simple cell-and-network pattern (plain line). Both models 
have the same average opacity and source function. The factor-2 
jump at 1/fi = 1 1.3 corresponds to the limb of the opaque solar 
disk; the reference radius used to compute ju corresponds to the 
middle of the emitting layer. The oscillations are artefacts of the 
averaging process. 



ficients k Vo . We then solve the radiative transfer equations for Aq 
along rays coming from infinity through the emitting layer to the 
observer. 

As the opacity is obtained by a simple integration of k V(j ns vl j, 
the average line-of-sight opacity to as a function of fi for the 
"paddle-wheel" pattern is the same as for a uniform layer with 
the same average fc Vo Hsvi,;- However, as seen in Fig. [7] still for the 
same average S and k Vo ns vl j, the effect of opacity (a decrease 
in intensity) is higher in the "paddle-wheel" case, in particular 
for intermediate values of 1 /fi. As a result, neglecting the cell- 
and-network pattern of the real TR leads to overestimating the 
opacity when using method A. 



5.4. Roughness and fine structure 

In order to expla in the high values of opacity (as derived from 
their method A), Dumo nt et al.l ([1983) introduce the concept of 
roughness of the TR: as the TR plasma is not perfectly vertically 
stratified (there is some horizontal variation), method A leads 
to an overestimated value of To- This could reconcile the values 
obtained following our application of methods A and B. 

We model the roughness of the transition region by incom- 
pressible vertical displacements of any given layer (at given op- 
tical depth) from its average vertical position, in the geometry 
shown in Fig. [8] The layer then forms an angle a with the hori- 
zontal and has still the same vertical thickness ds; the thickness 
along the LOS is ds cos a/ cos(6> + a), as can be deduced from 

Rg. m 

If we assume that 9 + a remains sufficiently small for 
the plane-parallel approximation to hold (and so that the LOS 
crosses one given layer only once), the opacity is 



to 



= fn t 



(s)K(T(s)) 



cos ads 
cos(9 + a) 



(24) 



The angle a is a random variable, with some given distribu- 
tion Pr(o'). We compute the average of to as a function of 9 and 



7 For high values of the Pr(o') width A, the correction factor /? cannot 
be computed for high values of l//j (high angles 9) because the values 
of a in the wings of Pr(a) fall in the range where \6 + a\ n/2: the 
plane-parallel approximation is not valid anymore. This explains the 
limited range of the /3(l/yu) curves in Fig. [9] 
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Figure 8. Geometry of a TR layer (plain contour), displaced 
from its average position (dashed contour) while retaining its 
original vertical thickness ds, and locally forming an angle a 
with the average layer. The line-of-sight (LOS) forms an angle 8 
to the vertical (normal to the average layer). 
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Figure 9. Multiplicative coefficient to fo due to roughness, for 
different roughness parameters A. 



ent methods from observations of the full Sun: center-to-limb 
variation of radiance, center-to-limb ratios of radiance and line 
width, and radiance ratio of the 93.3-94.4nm doublet. Assuming 
a spherically symmetric plane-parallel layer of constant source 
function, we find a S vi 93.3 nm opacity of the order of 0.05. The 
derived average electron density is of the order of 2.4 • 10 16 irT 3 . 

We have then used the line radiance (by an EM method) 
in order to get the RMS average electron density in the 
S vi 93.3 nm-emitting region: we obtain 2.0 ■ 10 15 irT 3 . This 
corresponds to a total pressure of 10~ 2 Pa, sligh t ly hi gher than 
the range of pressures found by Dum ont et al.l d 19831) (1.3 to 
6.3 • 10~ 3 Pa, as deduced from their Sec. 4.2), but lower than the 
value given in lMariskal (1 1993b (2 • 10~ 2 Pa). 

The average electron densities obtained from these methods 
(opacity on one hand, radiance on the other hand) are incompat- 
ible, as can be seen either from a direct comparison of the values 
of (n e ) and (n 2 ) for a given thickness As of a uniform emitting 
layer, or by computing the As that would reconcile the measure- 
ments of (n e ) As and (n 2 ) As. Furthermore, we have seen that 
the density obtained from the opacity method is also incompati- 
ble with standard DEMs of the Quiet Sun (see Sec. |4.2| > and with 
semi-empirical models of the temperature and density profiles in 
the TR (see Sec. 15X21 

We investigated several possible sources of biases in the de- 
termination of to: the approximation of a constant temperature 



in the S vi emitting layer, the anomalous behavior of the S vi ion, 
the chromospheric network pattern, and the roughness of the TR. 
Some of these could help explain partly the discrepancy between 
the average densities deduced from opacities and from radiances, 
but there is still a long way to go to fully understand this dis- 
crepancy and to reconcile the measurements. At this stage, we 
can only encourage colleagues to look for similar discrepancies 
in lines formed around log T — 5.3 (like C iv and O vi), Na-like 
and not Na-like, and to repeat similar S vi center-to-limb mea- 
surements. 

In Sec. 15.1.21 we have tried to combine opacity and radi- 
ance information to compute the gradient of temperature. This 
appeared to be impossible (if restricting ourselves to a realistic 
range of parameters) because of the above-mentioned incompat- 
ibility. We have estimated that a value tq of the S vi 93.3 nm 
opacity compatible with radiance measurements and with real- 
istic values of the temperature gradient would be in the range 
5 10~ 3 to 10~ 2 . 

In spite of the difficulties we met, we still think that the com- 
bination of opacity and radiance information should be a pow- 
erful tool for investigating the thermodynamic properties and 
the fine structure of the TR. For instance the excess opacity 
derived from observations and a plane-parallel model could be 
used to evaluate models of roughness and fine structure of the 
TR. Clearly, progress in modelling the radiative output of the 
complex structure of the TR needs to be done in order to achieve 
this. 
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Appendix A: About the filling factor 

In this paper we have defined the filling factor as 



/ = 



<» e ) 2 

<"e 2 > 



(A.l) 



while it is usually obtained, from solar observations (e.g. Judge 
2000; lKhmchuk & CargihfeOOll) . from 



EM 
AsiIq 



(A.2) 



where EM is the emission measure, As is the thickness of the 
plasma layer and «o is the electron density (usually determined 
from line ratios) in the non-void parts of the plasma layer. 

It may seem surprising that the EM is at the numerator of this 
second expression, while it provides an estimate for (« 2 ) which 
appears at the denominator of the first expression. However, we 
can show that these both expressions, despite looking very dif- 
ferent, give actually the same result for a given plasma. 

We take a plasma with a differential distribution £(n e , T) for 
the density and temperature: £(« e , T) dn e dT is the proportion 
of any given volume occupied by plasma at a density between n e 
and n e + dn e and a temperature between T and T + dT. 
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The contributions to the line radiance E and to the opacity at 
line center tq from a volume V with this plasma distribution are 



- = JJ n 2 e G(n e , T) f(n e , T) An, AT (A.3) 
^ = JJ « e ^(»e, r) f (He, 7") An e dr (A.4) 



with the notations of our article. 

The usual assumption (e.g. Judge 2000) is that G(T, n e ) "se- 
lects" a narrow range of temperatures around T = r max and does 
not depend on n e , i.e., G(n e , T) » G(T max ) 5(T - r max ). Similarly, 
we can consider that K(n e , T) m K(T max ) S(T - T max ). Then 



E 
V 

la 

v 



G(T m . dx ) J nl %(n e , T m . dx ) An e = G(T mRX ) (nl) T=Tmia 

(A.5) 

K(T mRX ) J n e £(n e ,T max ) An e = K(T mRX ) <n e >r=r, mx 



(A.6) 



The line ratio Rjj = Ej/Ej is, following [judge (2000]) and 
with the assumption G(n e , T) = G(n e , r max ) <5(T - r max ): 



_ _ JJ» e 2 G,(» e ,r)^(» e ,r) d» e dr 
~ £V " JJ«;k;/«e,r)£(«e,r) dn e dr 

J ne G '("e,r max )^(n e ,rmax) d« e 

J nlGj(n e ,T mRX )^(n e ,T mRX ) An e 

When homogeneity is assumed, i.e., £(« e , T) = 6{n e 
this becomes 



(A.7) 
(A.8) 

«o)l(n 



ngG/(no, r max )|(r max ) _ G,(«o, r max ) 
n^G /n , r max ) |(r max ) G/no, T max ) 



Sgy(fIo) (A.9) 



and inverting this function allows to recover «o from the ob- 
served value of Rjj. 

The fundamental point is that Rjj does not depend on the pro- 
portion / (the filling factor) of the volume actually occupied by 
the plasma: no is the density in the non-void region only. For ex- 
ample, for £f(n e , T) defined by fd(n e -«o) + (l ~f)8(n e ), the line 
ratio Rjj is gij(no) which is independent on /, while («e)r=r M 
determined from E/V would be fn^ and (n e )r=r max determined 
from Trj/V would be friQ. One can see in this case that / can 
(equivalently) either be recovered from 



<^>r=r ma = Cfrg) 
(no) 2 (no) 2 

(corresponding to Judge 2000) or from 



(nlh 



(/«o) 2 



n=T nm (fnf) 
(corresponding to our method). 



/ 



(A. 10) 



(A.ll) 
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